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Abstract 



The data collected at centre-of-mass energies of 188.6 GeV by ALEPH at LEP, 
corresponding to an integrated luminosity of 176.2 pb _1 , are analysed in a search for 
pair-produced charged Higgs bosons H 1 * 1 . Three analyses are employed to select the 
t~*~v t t~v t , cst~v t and esse final states. No evidence for a signal is found. Upper 
limits are set on the production cross section as a function of the branching frac- 
tion B(H + ^t + v t ) and of the mass mn±, assuming that the sum of the branching 
ratios is equal to one. In the framework of a two-Higgs-doublet model, charged Higgs 
bosons with masses below 65.4 GeV/ c 2 are excluded at 95% confidence level indepen- 
dently of the decay mode. 
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1 Introduction 



In the Standard Model, particle masses are generated via the Higgs mechanism implemented 
using one doublet of complex scalar fields. In this process one physical state remains 
in the spectrum, known as the Standard Model Higgs boson. The most important 
phenomenological consequence of an extended Higgs sector is the appearance of additional 
spin-0 states, both neutral and charged. For example, with the addition of one more doublet, 
five physical states remain after the spontaneous breaking of the SU(2)l x U(1)y symmetry: 
three neutral and a pair of charged Higgs bosons. Among the possible extensions of the 
Higgs sector, those obtained by adding more doublets are preferred because they naturally 
lead at tree level to M\y ~ Mz cos#w, a relation very well verified by experiment. 

The ALEPH data collected at centre-of-mass energies up to 184 GeV have been used 
in Refs. |], to search for pair production of charged Higgs bosons predicted in models 
with two Higgs doublets. The negative result of the search was translated into a lower limit 
on the H ± mass m^± of 59 GeV/c 2 at 95% confidence level. In this paper an update of 
the search based on the data collected at ^Js = 188.6 GeV (hereafter referred to as the 
189 GeV data) is presented. The theoretical framework and underlying assumptions are 
the same as detailed in Refs. |l|, ||. The H + is assumed to decay predominantly into cs or 
t + v t final states (and respective charge conjugates for the H~). Other decay modes are 
not considered and B(H + — »t + v t ) +B(H + — >cs) = 1 is assumed, but the analysis is equally 
sensitive to other hadronic decay modes. As a consequence, H + H~ pair production leads 
to three final states (esse, csT _ z/ T /csr + z/ r and t + v t t^v t ) for which separate searches are 
performed. 

This letter is organized as follows. After a brief description of the ALEPH detector in 
Section 2, the event selections are described in Section 3. The results and the conclusions 
are given in Sections 4 and 5. 

2 The ALEPH detector 

Only a brief description of the ALEPH subdetectors relevant for this analysis are given 
here. A more comprehensive description of the detector components is given in Ref. || and 
of the reconstruction algorithms in Ref. H. 

The trajectories of charged particles are measured with a silicon vertex detector, a 
cylindrical drift chamber, and a large time projection chamber (TPC). These are immersed 
in a 1.5 T axial field provided by a superconducting solenoidal coil. This system yields a 
resolution of Spx/pr = 6 x 10~ 4 />r©0.005 {pr in GeV/ c). Hereafter, charged particle tracks 
reconstructed with at least four hits in the TPC, and originating from within a cylinder of 
length 20 cm and radius 2 cm coaxial with the beam and centred at the nominal collision 
point, are referred to as good tracks. 

The electromagnetic calorimeter, placed between the tracking system and the coil, is 
a highly segmented sampling calorimeter which is used to identify electrons and photons 
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and to measure their energies. It has a total thickness of 22 radiation lengths at normal 
incidence and provides a relative energy resolution of 0.18 /y/E + 0.009 (E in GeV). The 
luminosity monitors extend the calorimetric coverage down to 34 mrad from the beam axis. 

Muons are identified by their penetration in the hadron calorimeter, a 1.2 m thick iron 
yoke instrumented with 23 layers of streamer tubes, together with two surrounding layers 
of muon chambers. The hadron calorimeter also provides a measurement of the energy of 
charged and neutral hadrons with a relative resolution of 85%/\/E~ (E in GeV). 

The calorimetry and tracking information are combined in an energy flow algorithm, 
classifying a set of energy flow "particles" as photons, neutral hadrons and charged particles. 
From these objects, jets are reconstructed with an energy resolution of (0.60\/~E + 0.60) x 
(1 + cos 2 9) where E in GeV and 9 are the jet energy and polar angle, respectively. 

3 Analysis 

To ensure good potential for discovery, independent of the branching fraction B(H + — >t + u t ), 
three selections are defined for the topologies t + v t t~v t , cst~v t /cst + v t (hereafter referred 
to as cst~v t ) and esse. The most relevant selection criteria are chosen to achieve the best 
expected confidence level for exclusion of a mass hypothesis of 70 GeV/ c 2 . Each selection is 
optimised individually with the most optimistic B(H + ^t + u t ) in each case, i.e., 0%, 100% 
and 50% for the esse, t + v t t~v t and cst~u t channels, respectively. 

3.1 Monte Carlo samples 

Fully simulated Monte Carlo event samples reconstructed with the same program as the 
data have been used for background estimates, design of selections and cut optimization. 
Samples of all background sources corresponding to at least 20 times the collected luminosity 
were generated. The most important background sources are e + e~ — > r + r~, qq, four- 
fermion processes and two-photon collisions, simulated with the KORALZ ||, PYTHIA ||, 
KORALW and PHOT02 § generators. 

The signal Monte Carlo events were generated using the HZHA || generator. Samples 
of at least 2000 signal events were simulated for each of the various final states for charged 
Higgs boson masses between 50 and 75 GeV/c 2 . 

3.2 The t + v t t~v t final state 

The final state produced by leptonic decays of both charged Higgs bosons consists of 
two acoplanar r's and missing energy carried away by the neutrinos. As this topology 
is identical to that expected from stau pair production with massless neutralinos, the 
"Large AM" selection described in Ref. (JTOj is used here to search for charged Higgs 
bosons in this channel. Efficiencies to select events from H + H~ — ► t + v t t~v t are of the 
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Table 1: Selection efficiencies e (in %) as a function of the charged Higgs boson mass mn±. 



m H ± (GeV/c 2 ) 


50 


55 


60 


65 


70 


75 


e (t + v t t~v t ) 


34 


35 


38 


35 


40 


41 


e (cst~9 t ) 


35 


37 


37 


35 


29 


20 


e (esse) 


48 


48 


49 


49 


48 


45 



order of 35%, as shown in Table [3] for a representative set of charged Higgs boson masses. 
The total expected background amounts to 15.5 events, mainly consisting of irreducible 
background from W + W~ — > t + v t t^v t . In the data, 20 events are selected, in agreement 
with the expectation. The systematic uncertainty on the number of expected signal events 
is estimated to be 3.0%, dominated by the effect of limited Monte Carlo statistics (2.7%) 
and uncertainties on the cross-section for charged Higgs production (1.0%). The systematic 
error on the background is estimated to be 8%. This is dominated by the effect of limited 
Monte Carlo statistics (4%), uncertainties on the cross section for W pair production (2%), 
and uncertainities on the cross section for two-photon processes (7%). The systematic error 
on the luminosity is estimated to be 0.5%. 

3.3 The cst~v t final state 



The mixed final state csr~z/ r is characterised by two jets originating from the hadronic 
decay of one of the charged Higgs bosons and a thin jet plus missing energy due to the 
neutrinos from the subsequent decay of the charged Higgs and of the r. 

As a first step in the analysis the thrust of the event is required to be less than 0.96 and 
the total number of good charged tracks greater than 7. To reduce background from two- 
photon processes and the contribution of beam related background which is not simulated, 
the energy deposited in a 12° cone around the beam axis is required to be less than 2.5% of 
the centre-of-mass energy. Background from e + e~ — > qq(7) events is reduced by demanding 
that the polar angle 6 m i SS of the missing momentum vector point away from the beam axis 
such that Icos^missl < 0.9. To reduce the background from W + W~ — ► qq'-ft^ where £ 
corresponds to an electron or muon, it is required that the events contain no identified 
lepton with a momentum greater than 10% of the centre-of-mass energy. 

At this point the events are clustered into three jets using the JADE algorithm [TT[. The 



Vent value where the transition from two to three jets occurs is required to be greater than 
0.001. The jet with the lowest charged track multiplicity is taken as the r jet candidate. 
If at least two jets have the same multiplicity the r jet candidate is taken to be the lowest 
momentum jet. The following quality cuts are then applied on the r jet. It is required that 
the charged multiplicity of the r jet be between one and three. The angle between the r 
jet and the closest quark jet candidate is required to be between 30° and 125°. The energy 
of the t jet boosted into the charged Higgs boson rest frame, defined as the frame recoiling 
against the hadronic system, is required to be less than 40 GeV. 

To further suppress backgrounds the following four variables are used: 
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• The angle 6^ between the two hadronic jets. 

• The total transverse momentum p\ ot divided by the total visible energy E V i S . 

• The production angle of the events, reconstructed from the sum of the momenta of 
the two quark jets. This production angle is then charge tagged with the tracks from 
the r jet to form the variable 0^ od . In the case of two tracks in the r jet the charge 
of the highest momentum track is used. 

• The x 2 P er number of degrees of freedom from a kinematic fit to the events using the 
constraints of energy and momentum conservation and the equality of the two masses 
produced in each side of the event. 



The four variables are linearly combined in a discriminating variable D, displayed in 
Fig. [l|a. The cut optimisation leads to D > 0.26. In the data collected at y/s = 189 GeV, 20 
events are selected, in agreement with the background expectation of 22.6. The efficiencies 
of selection for a range of masses are given in Table [I]. The mass of the hadronic system, 
and hence the mass of the charged Higgs candidates, is rescaled such that the energy of the 
two jets is equal to the beam energy in order to improve the resolution. The reconstructed 
masses of the candidate events are displayed in Fig. |l]b. The cutoff near 80 GeV/c 2 is due 
to the influence of the # q q/ variable. 

The systematic uncertainty on the number of expected signal events is estimated to 
be 4.1%. The main contributions are from Monte Carlo statistics (3.5%), calibration 
uncertainties (1.5%) and uncertainty on the cross section for charged Higgs production 
(1.0%). The systematic error on the background level is estimated to be 12%. The main 
contributions are from Monte Carlo statistics (3%), uncertainty on the cross section for the 
W + W~ process (2%) and from the statistics of data/Monte Carlo comparisons (10%). The 
systematic error on the luminosity is estimated to be 0.5%. 



3.4 The esse final state 

The hadronic decays of two charged Higgs bosons lead to a final state of four well 
separated jets, which can be combined into two dijets with equal masses. With respect to 
Ref. H the preselection and jet pairing method remain unchanged. 

The following five variables are used: 

• The x 2 P er degree of freedom from a five-constraint kinematic fit. The constraints in 
the fit are from conservation of energy and momentum and the equality of the two 
charged Higgs masses in the event. 

• The production polar angle 9 pro d, i.e., the angle between the charged Higgs boson 
momentum direction and the beam axis. 

• The difference between the largest and the smallest jet energies, E max — E min . 
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Figure 1: (a) The distribution of the discriminating variable D described in the text for 
the semi leptonic channel at the level of preselection, (b) The reconstructed masses of 
the charged Higgs candidates after the cut on the discriminating variable. In both plots 
the points are the data, the open histograms are the Standard Model background and the 
hatched areas represent the charged Higgs signal expectation with m H ± = 70 GeV/c 2 . The 
signal is normalised arbitrarily in both plots. 



The product of the minimum angle between any two jets, and the smallest jet energy, 



The QCD matrix element Aiaq' [12 



The variables are linearly combined into one discriminating variable D, shown in Fig. §a. 
Events are accepted if D > 4.4. Including in D the charm tag used in Ref. [Q] does not 
increase the discriminating power. 

The analysis selects 263 events from the data for masses between 50 and 80 GeV/c 2 , 
corresponding to a background of 294.2 events expected from Standard Model processes. 
The fitted-mass distribution of the selected candidates can be seen in Fig. ||b. Efficiencies 
are of the order of 50% as shown in Table [l|. 

The systematic error on the number of events expected is estimated to be 3.0%. The 
main contributions are from Monte Carlo statistics (2.4%), statistics of data/Monte Carlo 
comparison (1.2%), and knowledge of signal cross sections (1.0%). The systematic error on 
the background level is estimated to be 3%. The main contribution is from knowledge of 
the W + W~ cross section (2%). The contribution of the luminosity is estimated to be 0.5%. 
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3 4 5 6 50 60 70 80 90 100 

D m™ ± (GeV/c 2 ) 

Figure 2: (a) The distribution of the discriminating variable D described in the text for 
the esse channel at the level of preselection, (b) The distribution of the Higgs candidate 
masses after the cut on the discriminating variable. In both plots the points are the data, 
the open histograms are the Standard Model background and the hatched areas represent 
the Higgs signal expectation with m H ± = 70 GeV/c 2 . The signal is normalised arbitrarily 
in both plots. 



4 Results 

The numbers of candidates observed in the data collected at a centre-of-mass energy of 
189 GeV are consistent with those expected from Standard Model processes for each of the 
three channels. Since, in addition, the mass distributions in the esse and csr~z/ T channels 
do not show any significant accumulation outside the W + W~ region (Figs. [L] and 0), the 
results of the three selections described in this letter are combined with those obtained using 
y/s = 172-183 GeV data to set a 95% confidence level upper limit on the cross section for 
pair production of charged Higgs bosons. 

In setting the limits several new features with respect to Refs. fl], || are to be noted. 
Full background subtraction is performed according to Ref. |13|]. The likelihood ratio test 



statistic is used. The confidence levels are calculated using the semi-analytical approach 
described in Ref. ||14]| . Systematic errors are conservatively taken into account by reducing 



the efficiencies and subtracted backgrounds by one standard deviation. The reconstructed 
mass of the charged Higgs boson is used as discriminating variable for the esse and csr~z/ T 
channels. 

The upper limit on the H + H~ production cross section at 188.6 GeV as a function of 
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Figure 3: Upper limits at 95% confidence level on the H + H~ production cross section at 
s/s = 188.6 GeV for three values of B(H + ^t + v t ). The charged Higgs boson production 
cross section is shown solid curve. 



m H i is shown in Fig. |3| for three values of B(H + — ►r + z/ r ). The results from lower centre- 
of-mass energies have been scaled to -^/i = 188.6 GeV according to the dependence of the 
cross section on the centre-of-mass energy. 

In two-Higgs-doublet models the production cross section for H + H~ depends, at lowest 
order, only on mi±. The expected cross section at 188.6 GeV, corrected for initial state 
radiation, is shown in Fig. [| Upper limits on the production cross section therefore translate 
into excluded domains for mn± . 

The result of the combination of the three analyses is displayed in Fig. |j. Charged 
Higgs bosons with masses less than 65.4 GeV/c 2 are excluded at 95% confidence level 
independently of B(H + ^t + u t ). The corresponding expected exclusion is 69.1 GeV/c 2 . 
For the values B(H + -^T + v r ) = 0, 0.5 and 1, 95% CL. lower limits on m H ± are set at 69.9, 
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65.4 and 76.3 GeV/ c 2 , respectively. 



5 Conclusions 

The search for pair-produced charged Higgs bosons in the three final states t + u t t~v t , 
csr~z/ T and esse has been updated using 176.2 pb -1 of data collected at y/s = 188.6 GeV. 
No evidence of charged Higgs boson production has been found and upper limits have been 
set on the production cross section as a function of B(H + — ^ t + v t ) and of m H ±. Within 
the framework of two-Higgs-doublet models these results exclude at 95% confidence level 
charged Higgs bosons with masses below 65.4 GeV/c 2 independently of B(H + ^t + v t ) and 
assuming B(H + ^r + z/ T ) + B(H + ^cs) = 1. Similar results have been reported by L3 |0| . 
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Figure 4: Limit at 95% C.L. on the mass of charged Higgs bosons as a function of B(H + -^t + v t ). 
Shown are the expected (dash-dotted) and observed (solid) exclusion curves for the combination 
of the three analyses, and the full 172-189 GeV data set. The shaded area is excluded at 95% 
C.L.. The observed limits from individual channels are indicated by a thin solid line. 
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